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Abstract 

Understanding random lasing is a formidable theoretical challenge. Unlike conventional lasers, 
random lasers have no resonator to trap light, they are highly multimode with potentially strong 
modal interactions and they are based on disordered gain media, where photons undergo random 
multiple scattering [IH2] • Interference effects notoriously modify the propagation of waves in such 
random media, but their fate in the presence of nonlinearity and interactions is poorly understood. 
Here, we present a semiclassical theory for multimode random lasing in the strongly scattering 
regime. We show that Anderson localization [HIS], a wave-interference effect, is not affected by 
the presence of nonlinearities. To the contrary, its presence suppresses interactions between simul- 
taneously lasing modes. Using a recently constructed theory for complex multimode lasers |6], we 
show analytically how Anderson localization justifies a noninteracting, single-pole approximation. 
Consequently, lasing modes in a strongly scattering random laser are given by long-lived, Ander- 
son localized modes of the passive cavity, whose frequency and wave profile does not vary with 
pumping, even in the multi-mode regime when mode overlap spatially. 



1 



Wave interference effects influence a broad range of physical processes in wave propaga- 
tion through random media [3 |H], ranging from the time-evolution of ocean waves and the 
transmission of light through interstellar clouds, down to the dynamics of individual elec- 
trons in nano-electronic circuits and the transmission of photons through imperfect optical 
fibers and other waveguides. They generate a rich variety of wave-interference phenomena 
such as Anderson localization [H [5] and coherent backscattering [914TT] , as well as effects of 
higher-order coherence such as quantum noise and particle bunching or anti-bunching [T2] . 
For linear wave equations and noninteracting fields, these effects are well understood [H], 
however little is known of the influence that nonlinearities and interactions have on coherent 
wave interference effects. In random lasers, confinement of light within an amplifying media 
is achieved via multiple scattering P-[3] and one thus expects important wave-interference 
effects to be present. The mechanism behind the observed spectral narrowing has been and 
still is a matter of debate, with most theoretical proposals constructing lasing modes from 
modes of the passive cavity trapped by scattering, be they Anderson localized modes in 
the strong scattering regime [T3] or anomalously trapped modes in the weak scattering 
regime [15]. Numerical calculations suggest that lasing in weakly scattering random lasers 
can be supported by extended diffusive modes flQ\ ITf]. Recent experiments have observed 
delocalized lasing modes in weakly scattering random lasers [HI [19], while lasing modes 
were significantly less extended for stronger scattering [20], with some evidence for mode- 
dependent localization lengths [21J. The Anderson localized regime has been probed in few 
experiments, where the spatial extent of lasing modes could not be directly accessed [22] ■ 
Because random lasers mostly operate in the multiple scattering, multi-mode regime, they 
are perfect testbeds to investigate effects of wave coherence in the presence of nonlinearities 
and interactions. With the exception of Ref. [17], theoretical and numerical investigations of 
the localization properties of lasing modes rely on calculations of the total electromagnetic 
field inside the amplifying medium. There is still very little understanding of the local- 
ization properties of individual lasing modes, how they are influenced by the infinite-order 
nonlinearity arising from spatial-hole burning and by the resulting modal interactions. Of 
particular interest is to find out if, above lasing threshold, modal interactions change the 
localization length in an Anderson localized passive cavity. The influence of modal inter- 
actions usually becomes more and more pronounced as more modes start to lase, thus it is 
expectable that lasing modes differ more and more from passive cavity modes as the pump 
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strength increases. 

We focus on the interplay between Anderson locahzation and nonhnearities in a multi- 
mode random laser, which we take as a photon-scattering gain medium region surrounded 
by vacuum. Below threshold, the system is linear and we assume the gain medium to be 
strongly scattering, so that the modes of the passive cavity are exponentially localized over a 
localization length ^ much smaller than the cavity's linear size, ^ ^ L. This simultaneously 
means that modes are long-lived, with a mean lifetime ^ / c) exp[L / ^] and low lasing 
thresholds when they are localized deep enough inside the cavity [23] . Our goal is to under- 
stand how lasing modes individually differ from passive cavity modes. To that end we adopt 
the method of Ref. [6| (see Supplementary Information), where the electric field is expanded 
as E{x,t) = exp[— ifi^t] over individual lasing modes \E'^. The latter are in their 

turn expanded over so-called constant flux (CF) modes, which are eigenmodes of the passive 
cavity satisfying the boundary condition that they couple only to undamped outgoing waves 
(i.e. with real wavenumber k^) outside the cavity. The resulting modal resolution of E{x) 
allows to directly compare individual lasing modes with the CF modes. Our main finding 
is that, quite surprisingly, simultaneously lasing modes in the Anderson localized regime do 
not interact with each other. To a very high precision, the frequency and wavefunction of 
a lasing mode are the same as that of a CF state, regardless of the pump strength, even 
in the multi-mode regime and as new modes start to lase. Nonhnearities only decide which 
mode starts to lase at what pump threshold. Our analysis of the lasing theory of Ref. [H], 
augmented by wavefunction correlations in the Anderson localized regime [25, provides a 
straightforward analytical understanding of the suppression of modal interactions: a single- 
pole approximation is rigorously justified, even in the multi-mode regime, because of the 
exponentially small broadening of localized modes and their exponentially small coupling 
to only an algebraic number of CF modes. The absence of modal interaction we predict is 
totally unexpected: a similar system, with the same gain medium and cavity shape and size, 
but without or with little scattering standardly exhibits strong mode competition, with in 
particular shifts in mode frequencies and the disappearance (and possible re-appearance) of 
certain modes as the pump strength increases [T^. Modal interactions are turned on once 
CF modes are broadened, which can occur either at weaker disorder when their localization 
length increases, or when their lifetime is limited by photon absorption. The linear nonin- 
teracting behavior of lasing modes we predict is the trademark of Anderson localization in 
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random lasers, and we propose to monitor lasing frequencies as a function of pump strength 
to detect Anderson localization in experiments on random lasers. 

As starting point we take the Maxwell-Bloch equations (MBE), a set of three coupled 
nonlinear equations for the electric field, the polarization and the population inversion in 
the gain medium [25] (see Supplementary Information). For simplicity we consider either a 
one-dimensional system or a two-dimensional system with electromagnetic field transverse 
electric with polarization perpendicular to the gain medium, or transverse magnetic. In those 
instances, it is sufficient to consider scalar fields. To find steady-state solutions, Ref. [H] first 
writes both the polarization and the electric field in a multiperiodic expansion. Further 
imposing the rotating wave and stationary inversion approximations, one obtains a set of 
self-consistent equations 

DoTa^ = DoAra^, (1) 



a 



for the vector a^ = (a^, a^, ) of components of the expansion (we use dimensionless units 
for ^ and Dq, as described in the Supplementary Information) 

vl/^(x,fc,) = ^a;r7A™(x,fc^), (2) 

m=l 

of the lasing modes \E'^(x, /c^) over elements of a truncated basis of CF modes 
{iprn{^: ^M)}m=i- When modes are lasing, M > N, and in the Anderson localized regime, 
we will see that the structure of the threshold matrix T is such that = 1 for each lasing 
state. In Eq. ([T]), Dq gives the overall pump strength, which we define as the population 
inversion density of the active medium. The threshold matrix T = AT is the product of a 
diagonal and a non-diagonal matrix. 



Am„ - 5mn ,n I'n , _ , . , , (3a) 



Tmn = I d'^x— ^^^^Vm(x, A;^)^„(x,fc^) , (3b) 

with the pump profile -F(x), the speed of light c, the atomic frequency ujq of the two-level 
active medium and the polarization relaxation rate 7_l. The integral runs over the volume 
C (length in dimension ci = 1, area in ci = 2) of the gain medium. Nonlinearities arise from 
the hole burning denominator 

Mx) = 5^r(M|M/.(x)p, (4) 



where the sum runs over all lasing modes and T{k) = 'yj_/[{ck — wq)^ + 71]- The lasing 
modes are self-consistently determined by the nonlinear eigenvalue problem of Eq. ([T]). So- 
lutions exist only for discrete real values of k^, giving the lasing frequencies = ck^. The 
wavenumber of the CF modes km^i = km^ + ii^mfi, with fcm^, t^^n^l ^ Ti, are complex because 
of the openness of the cavity, k"^ = CTm^ with the lifetime Tm^. 

The presence of modal interactions is reflected in the structure of the lasing modes, 
e.g. the number of CF states significantly contributing to a single lasing mode in the 
expansion of Eq. (|2]). If the expansion is dominated by a single component, the lasing mode 
resembles a noninteracting CF mode, however if the expansion runs over many nonvanishing 
components, the lasing mode is in general very different from noninteracting CF states, 
signaling the onset of modal interactions. We next proceed to show that ~ 1 in the 
strongly localized regime, by analyzing the structure of the threshold matrix T = AT. 
In the Supplementary Information, we argue that matrix elements %nn are non-zero only 
when CF states ipm and ipn are centered within one localization length of one another. 
Furthermore, the nonvanishing matrix elements fluctuate pseudo-randomly as a function 
of the indices (m, n). In the Anderson localized regime, most CF modes are long-lived, 
^^m^l - ^~^exp[-L/^] and consequently \Kmm\ - (fc^^7±/2) exp[L/^] /{ck^ - (^o + i7±) is 
exponentially large at k^ = km^i- One then expects that lasing frequencies are very close to 
those of certain CF modes, k^ ~ k^ofj., and that the expansion ^ for the corresponding 
lasing mode is dominated by a single component aj^o. We show perturbatively that this is 
indeed the case as long as exp[— L/^] <^ 1. Keeping only a single component in Eq. ([T]) gives 
the single-pole approximation, {a^Y^^ = for n 7^ mo, and 

(a7)(°) = DoT™,^Ja7)W. (5) 

At this level, one has Do = Tmlmo = i^momoTmomo)'^ oc exp[-L/^] < 1. Leading-order 
corrections are obtained with the approximation 

From Eqs. ([s]) and ([6]), we conclude that corrections to the single-pole approximation are 
negligible when 

^ _ SnT^mo ^nn%non%imo , . 

^= 7 [') 

^^momo 'mo mo 
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with Kmn and Tmn defined in Eqs.([3]). We next analyze the conditions under which Eq. ([T]) 
is satisfied. 

In order for Tmon%imo not to vanish, the CF states mo and n must have their support 
within one localization length. At energy hck^ they must thus differ in energy by at least 
the CF mode spacing within one localization volume pJJ, in our case = hc/bdk'^~^^'^, 
hi = 2 and 62 = 27r. We have the upper bound |A„„| < h'^c^kf^/ {2hck^/\^ + A|) x \^^/{ck^ — 
^0 + ^7±)l fo^^ ^ 7^ ^0 in Eq. ([t]), and thus \Ann/ J^momo 

2hd{k^^Y]. Next, averages TmonTnmo of product of matrix elements of T can be evaluated, 
under the assumption that CF states are sufficiently close to localized modes of a disordered 
closed system. This assumption is satisfied when ^ L. We obtain (see Supplementary 
Information) < /(A;^^) exp[— L/^], with f{x) an algebraic function of x. Furthermore, 
it can also be shown that fiuctuations of S are also exponentially small, oc exp[— L/^]. This 
gives an exponentially small upper bound for corrections to the single-pole approximation, 
both for d = 1 = i) and for d = 2 when ^ is not exponentially larger than the wavelength 
2'K/kfj_. We conclude that the condition expressed in Eq. ([T]) is satisfied and that a single- 
pole approximation is justified in the Anderson localized regime. Accordingly, we propose 
that the criteria for the absence of interactions in random lasers are ^ L in d = 1 and 
^ ^ L, kp^^ < L/^ m. d = 2. While we expect that the c? = 3 criterion is the same as in two 
dimensions, we recall that an analysis of the MBE for vector fields is required in this case. 

Corrections to the single-pole approximation become non-negligible when the mode life- 
time becomes shorter. We briefiy discuss two cases when this happens. First, in the Anderson 
localized regime, CF modes have an exponentially long lifetime, but their broadening can 
be limited by photon absorption. This can be included in the above estimate via the substi- 
tution Kmf, ^ ^"^exp[-L/^] Kra^, ^ {cT^hs)~^ ■ This leads to |A„„/AmomJabs < h/A^Tsi,s 

and corrections to the single-pole approximations become of the order S ~ fik^^)/ (A^^CTabs)- 
We conclude that photon absorption turns on modal interactions in the localized regime, 
when it reduces the mode lifetime such that Tabs ^ fik^^)/k^c. Second, mode lifetimes are 
shorter and wavefunction correlations are different in the d = 2 diffusive, delocalized regime, 
where all CF modes overlap. From Ref. [26], we estimate (see Supplementary Information) 
S ~ hf{k^£, k^L) / AiTmo^, with the level spacing = hc/2TTk^L'^. Prelocalized states 
exist in such systems [IS], whose lifetimes are long enough, Tprdoc ^ ^/^l that deviations 
from the single-pole approximation are still negligible if f{k^i, k^L) <^ AiTpreioc/^- Ref. [T7j 
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on the other hand suggests that lasing in diffusive random lasers can be supported by stan- 
dard generic diffusive modes. The hfetime of such modes can be estimated by the classical 
diffusion time td = L"^ /cd. through the system, and the dimensionless measure of corrections 
to the single-pole approximation becomes H ~ k^d x f{k^l, k^L) > 1. This agrees with the 
observed strong modal interactions reported in Ref. [TTj for a diffusive laser with ^ 1. 
We caution however that our estimates all rely on the assumption that CF modes have 
the same wavefunction correlations as modes of closed, nonabsorbing Anderson localized 
systems. This is no longer true in the weakly scattering, delocalized regime. 

We numerically check our theoretical predictions. From the above analysis, dimensional- 
ity plays only a limited role in the Anderson localized regime, and for numerical convenience, 
we take a Id edge-emitting, disordered microcavity laser. The cavity extends from x = —L/2, 
where we impose reflecting boundary conditions to x = L/2 where it is open. We fix its 
length, L = 2 fim and the resonant frequency Uq = ck^ with koL = 50. The transversal 
and longitudinal relaxation rates are 7± = 8 x IO'^Wq and 7|| = 2 x 10~^ci;o. Disorder is 
introduced in the cavity by a spatially varying refractive index randomly distributed in the 
interval n{x) G [1.3, 6.7] with a spatial correlation x = 8 t^^- This rather large variation has 
been chosen to ensure that with the numerically reachable values of L and ko, we are in the 
Anderson localized regime. The localization length of cavity modes is numerically estimated 
with the inverse participation ratio 

2 



dxl'^ix) 



dx\^{x)\\ (8) 



which estimates the spatial extent of a wavefunction regardless of its overall normalization. 
For an exponentially localized mode = exp[i9{x)] exp[—\x\/2C,]/y/2^, one has = 4^. 

In Fig. [T] we show a typical set of data for lasing modes in the strongly disordered, 
localized regime, indicating the mode's position, spatial extent, lasing frequency and output 
power. We estimate ^ ~ 0.03 — 0.1 /im from Fig. [l^, so that ^/i^ — 0.015 — 0.05 ^ 1 and 
we are in the Anderson localized regime. As usual, localization lengths vary from one mode 
to another. Lasing modes have long lifetimes, and are localized well inside the cavity, a 
distance larger than ^ away from the emitting cavity boundary. We see that as the pump 
varies by a factor of a thousand and up to twelve modes lase, lasing modes change neither 
their position (monitored by their center of mass. Fig. [i]d), nor their frequency (Fig. [T]i), 
nor their spatial extent (Fig. flk). This is even more remarkable, given that several modes 
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Figure 1: Evolution of lasing modes as a function of pump strength Dq for a non- absorbing, 
Anderson localized, half-open Id microcavity laser. The index of refraction is randomly distributed, 
n(x) € [1.3,6.7], with a correlation length x = 8 nm. In each panel, lines appear when the 
corresponding mode starts to lase. Panel a): Inverse participation ratio defined in Eq. ([s]), giving 
a measure of the spatial extent of the mode. Panel b): Position of the center of mass xq = 
J dx / dx of the mode. Panel c): Modal output power. The output power does 

not saturate because we did not limit the population inversion density in our numerics in order 
to investigate multi-mode lasing. Panel d): Lasing frequency in units of the population inversion 
relaxation rate 7_l. 

strongly overlap either spatially (for instance the red and brown modes) or in frequency. 
Only the modal emission power increases with pump strength, with a linear dependence on 
-Do that is reached quickly after threshold for each lasing mode. This is what is expected for 
noninteracting modes, given that Dq gives the inversion density of the active medium and is 
not bound in our approach. Incorporating a finite maximal density in our model would lead 
to mode saturation, but would also limit the number of lasing modes and thus constrain our 
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X [|j,m] X [^m] 

Figure 2: Spatial mode profile of lasing modes for a cavity with the same parameters as 

in Fig. [l| Left panel: at pump strength Dq = 0.04, there are seven lasing modes. Straight lines are 
exponential fits giving the indicated localization lengths Fits around x = —0.2fim (not shown) 
give ^ ~ O.OG^m. Right panel: at pump strength Dq = 0.4, three new modes have started to lase 
(dashed lines), however, lasing modes have not changed their overall profile, but have only been 
multiplied by a mode-dependent factor. 

investigations of mode interactions. 

We investigate more quantitatively the evolution with pumping of the spatial mode struc- 
ture. We show in Fig. [2] the spatial profile of lasing modes at two different pump strengths, 
both in the multi-mode regime. At moderate pump strength, Dq = 0.04 (left panel) there are 
seven lasing modes which all exhibit exponential Anderson localization. From exponential 
fittings, we have extracted localization lengths in the range ^ G [0.05, 0.3]/xm, in good agree- 
ment with the values for the inverse participation ratio shown in Fig. [1^. Upon increasing 
the pump strength to Dq = 0.4 (right panel), three additional modes are lasing (they are 
indicated by dashed lines), however the shape of the lasing modes has not changed at all. 
Enhanced pumping only increased their overall magnitude by a mode-dependent factor. 

Fig. |3] next shows in color scale the normalized spatial profile |\E'(x)p/ da;|\E'(x)p of 
the second lasing mode, corresponding to the red curves in Figs. [T] and |2} It is clearly seen 
that over almost four orders of magnitude for Dq, and despite nine mode thresholds being 
crossed, the mode's spatial structure remains exactly the same. Particularly interesting is 
the fact that around Dq = 0.3, a new mode starts to lase (brown curves in Figs. [l]and|2]), 
spatially right on top of the mode plotted on the left panel. Despite their strong overlap, 
these two modes do not interact, which we attribute to strong localization. 
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Figure 3: Color scale plots of the normalized spatial mode profile dx|^'(x)p for the 

second (left) and third (right) lasing mode in Fig. [T] corresponding to the red and green curves 
there, respectively. The arrows show modes thresholds, placed at the centers of mass of the modes, 
with the same color coding as in Figs. [T] and [2j It is clearly seen that the spatial structure of the 
modes is insensitive to the pump strength, even for the second lasing mode (left, corresponding 
to the red curves in Fig. [T|, which strongly overlaps with the brown mode starting to lase at 
Do ~ 0.03. 

Having confirmed our prediction that modal interactions are frozen in the Anderson 
localized regime, we finally illustrate how they are turned on when the lifetime of the passive 
cavity modes is reduced. To achieve this, we first add some absorption in the cavity via 
a finite imaginary part of the index of refraction. We show in Fig. and b data for the 
same cavity as in Fig. [T], where however the index of refraction acquires a finite imaginary 
part, lm[n] = 0.15. We see that modal interactions are turned on as inverse participation 
ratios and modal emission power vary significantly with the pump strength. We also found, 
but do not show, that the lasing frequencies and the centers of mass of the modes vary 
with the pump strength. We finally corroborate our prediction that modal interactions 
are turned on in the delocalized regime by simulations of cleaner cavities in Fig. ^ and d. 
We first observe that less modes lase, which we attribute to the presence of stronger mode 
interactions. The presence of strong interactions is further confirmed by the appearance and 
disappearance of a lasing mode (blue curves) as the pump increases. Such switching on and 
off of lasing modes with increasing pump strength have been reported in two-dimensional 
weakly disordered cavities [17] , where they were due to frequency repulsion. This is not the 
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Figure 4: Spatial extent, as measured by the inverse participation ratio (left), and modal output 
power (right) of lasing modes, as a function of pump strength Dq, for a non-absorbing, Ander- 
son localized, half-open Id microcavity laser. Panels a) and b) are for an absorbing cavity with 
Re[n(x)] G [1.3,6.7] and Im[n(x)] = 0.15. Panels c) and d) are for a non-absorbing cavity with 
weaker disorder, n(x) € [1.3,4.2] and extended modes. All other parameters are the same as in 
Fig.[T} 

mechanism at work here, as we found that the blue mode in Fig. |4|:; and d is separated by 
about 0.57_L from the (yellow) lasing mode with closest frequency. 

The interplay between nonlinear interactions and quantum/wave coherent effects in sys- 
tems with complex scattering is not well understood. Here, we have shown how strong 
localization effects freeze mode interactions in multi-mode lasers, even between modes that 
have significant spatial or frequency overlap. Our theory shows that Anderson localization 
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justifies a single-pole approximation where modes do not interact, and nonlinearities due to 
spatial-hole burning only determine which modes lase and with what threshold. Our results 
suggest a novel approach to investigate Anderson localization in strongly scattering random 
lasers by observing how much lasing frequencies vary with pump strength, which does not 
require to spatially resolve the lasing field. 
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Supplementary information 



I. STEADY-STATE AB INITIO LASER THEORY 

We give a brief overview of the theory developed by Tiireci et al. [Si] for steady-state 
multi-mode lasing. Its starting point is the semiclassical theory of lasing, where the electro- 
magnetic field is described by the classical Maxwell's equations, and is coupled to the quan- 
tum mechanically described polarization of the gain medium. We take the latter to be 
a collection of two-level atoms. The result is a system of three nonlinear coupled partial 
differential equations [S2] 

d'^^E^ = n(x)-2 [V'£;+-47r5fP+] , (Sla) 
diP+ = -{iujQ + -i^)P+ -ig^h-^E+D , (Sib) 
dtD = -i^\(D - D)^2ih-^[E+{P+Y - P+{E+y]. (Sic) 

We restrict ourselves to one- dimensional or two-dimensional transverse polarization and ac- 



cordingly consider only scalar field. In Eqs. (SI), the population inversion D(x, t) generates, 



in the presence of an electric field £'(x, t), a nonlinear polarization P(x, t) of the atomic 
medium, which in its turn, is coupled to D through E. Here, fiwo gives the spacing of the 
atomic levels, 7|| their relaxation rate and 7_l the relaxation rate for the polarization, -D(x) 
is the pump strength and n(x) is the cavity's index of refraction, related to the dielectric 
function by e = ra^. Finally, the dipole matrix element g couples the electric field with the 
population inversion, and we used the standard procedure of writing the electric field and 
polarization in terms of their positive and negative frequency components, P = P^ + P~ , 
E = E^ + E~ , neglecting the coupling of negative to positive components in the spirit of 
the rotating wave approximation. 



The set of Eqs. (SI), augmented by appropriate boundary conditions, describe nonlinear 
lasing. They cannot be solved exactly, and the standard procedure is to time-evolve them 
numerically long enough to find their steady-state solution. SALT is an alternative which 



determines steady-state solutions of Eqs. (SI ) with a much lesser computational price to pay. 



This of course assumes that steady-state solutions exist. We take them as multiperiodic, i.e. 
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expand the fields and 

E+{x,t) = ^^^(x)expH^]^^] (S2a) 

P+{x, t) = J2 W exp[-iQ^t] , (S2b) 

over lasing modes, = 1,2, ...N. SALT then uses the stationary inversion approximation, 
dtD = 0, to determine the components of the expansion as solutions of a self-consistent 
equation, 

(v^ + "^(x)^-a *.(x) - -*is;3^««r^*..(-) . (S3) 

where we parametrized the pump using a volume unit Vq = 5'^/eo^7± as -D(x) = DqF{-x)/Vo, 
splitting it into an overall dimensionless scale Dq and the spatial modulation -F(x), with 
maXx-F(x) = 1. The non-linearity source (the total weighted electric field intensity, or the 
hole burning denominator) is 

with the Lorentzian weight factors r(A;) = 7^/[(cA; — wq)^ + 7^], and electric intensity units 
^0 = /i^7lW/2^- 

Lasing theories usually construct lasing modes from resonances or quasiresonances of the 
corresponding (nonlasing) cold cavity. The latter are however non-normalizable outside the 
cavity, so that neither the modal power nor the relative weight of individual lasing modes 
can be extracted. SALT instead expands \l/^(x, A;^) = -E'o Xlm '^/TV'mlx, fc^), over a basis of 
so-called constant-flux ■i/'^ states which are defined by 

[V^ + n\y.)kl^^{k;)\ ^^(x, A:^) = 0, x G C , (S5a) 
[V' + n^fc^] ^^(x, A;^) = 0, x ^ C , (S5b) 

with the domain C defining the amplifying medium, rio the index of refraction outside and a 
requirement of containing only outgoing waves at infinity. The latter requirement, together 



with Eq. (S5b) can be recast as a condition on the constant flux state at the boundary of a 
d-sphere enclosing the resonator [d^ is the normal derivative at the surface of that sphere, 
d is the system's dimension) 
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As a consequence of this non-trivial boundary condition, the CF basis is not orthogonal in 
the usual sense. Instead, the following bi-orthogonality relation holds, 

j d'^Xn'^{x)i)*^{x,k^)lpn{^,k^) = 5mn, (S7) 

between the CF states and their duals ipj^{'K,kn). The latter are defined as solutions of 



complex conjugates of Eqs. (S5a) and (S6), which gives '^/'^(x, k^) = ^. 



Using the ansatz for the electric field intensity in Eq. (S3), one obtains a set of 



self-consistent equations for the vector of coefficients = (a^,a^, ...), 

a^ = DoTa^ = Z}oAra^. (S8) 

Here T is the threshold matrix, which depends on the set of vectors {a^}^^^, and which was 
split into the following two parts 

-k^ 



J^mn — ^mn ,n , , , . , , (S9a) 



Tmn = [ d'^x— ^^^y— ?/'^(x, fc^)?/'„(x, fc^) . (S9b) 
To calculate the modal power of the lasing state fi, we use a formula derived in Ref . |S3j . 



With our notation, it reads in c? = 1 dimension 



P, = 2e,Eluj,Sl ^ \^,{x)/E,\' (^T{k,)Doj^^^ - lm[n\x)]^ , (SIO) 

where eo is the vacuum permeability, and S is the cavity transversal cross-section. In the 
main text, we plot the modal power in dimensionless units, P^/{2eoEQUf^Sl), with / = 53 
nm. 



II. NUMERICAL APPROACH 



Solving Eq. (S8) self-consistently is, in general, highly non-trivial. The unknown vectors 



a^ enter the threshold matrix elements through the hole burning denominator h{x.) [see 



Eq. (S9b)] effectively leading to an infinite-order non-linearity. In addition, given a pump 
strength Dq, there is no simple way to find out how many modes are lasing, which is a 
prerequisite information for any numerical algorithm searching for a self-consistent solution 



to Eq. (S8). To construct the solution, we use an iterative scheme. We first determine 
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D=0 



Figure SI: Sketch of the microcavity laser we consider, where a Id cavity is closed to the left and 
open, hence emitting, to the right (dashed line). The pump area is shadowed. The spatial profile 
of a lasing mode (red) for a uniform pump and a spatially constant refractive index is plotted for 
illustrative purpose. 

the first lasing mode. Its pump threshold is entirely determined by the passive cavity 
properties — it is the mode of the largest lifetime in the absence of pumping. From here, 



the pump is increased by small steps. At each step, Eq. (S8) is iterated until convergence 



is achieved. Finally, the appearance of a real eigenvalue of the threshold matrix surpassing 
a certain level signals that a new mode starts to lase. This method has been worked out 
in a series of works by Tiireci, et al..|Slj and has been shown to predict lasing modes 
and threshold in complete agreement with the computationally much more expensive time- 



evolution of the Maxwell-Bloch equations, Eq. (SI ).|S4I 



The microcavity we model in our numerics is depicted in Fig. [STj It is closed on the left 
end, open on the right end, and pumped uniformly, F{x) = 1. In our numerics, we vary 
the degree of disorder in the refractive index, both in its strength and correlation length x- 
The disorder sets the CF states localization length ^. Typical values in our simulations are 
~ 100 — 1000 nm, A ~ 50 — 110 nm, and x = 8 so that ^ ^ A ^ x- Here A = 2t{ /uqUq 
is the wavelength of a freely propagating light wave in a material with a constant index of 
refraction, n{x) = no. 

Fig. |S2] illustrates structural changes of lasing modes in the interacting case. It should be 
compared with Fig. 3 of the main text, where a lasing mode is stable even in the exponentially 
small tails, despite the presence of many other lasing modes. Here we added finite imaginary 
part to the refractive index, to turn on modal interactions. They result in clearly visible 



profile changes of the state plotted in Fig. S2i: after the mode starts to lase, it tends to 
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1 10 100 1 10 100 

DO DO 

Figure S2: Color plots of the spatial mode profile ln|^'(x)p of the 2"^^ (left) and 13*^ (right) 
lasing mode for Id edge-emitting, disordered microcavity laser of Figs. 4a-b of the main text 
(corresponding to the red and maroon curves there). Colored arrows indicate the center of mass 
positions and pump thresholds of all other lasing modes, with the same color coding as in Fig. 4a-b. 

delocalize (its weight grows in the right half of the cavity). This growth is stopped as soon 
as other modes are turned on in that part of the space and the plotted mode mainly tends to 
localize more with enlarging the pump. The opposite behavior towards a weaker localization. 



is visible on Fig. p2p for another mode. These two examples suggest that there is no generic 
trend, and that both localization and delocalization may occur as the consequence of modal 
interactions. 



III. STRUCTURE OF THE THRESHOLD MATRIX 



The matrix elements of Tmn given in Eq. (S9b) are given by the pump profile, the spa- 



tial hole burning denominator and CF wavefunctions. We analyze these matrix elements, 
neglecting correlations between /i(x) and ip^i'^j fc^t)V'n(X) ^^t) = '^m(x, fc^)?/'„(x, /c^). We eval- 
uate averages {Tmn) and correlations {TmnTpq) of these matrix elements over an ensemble 
of cavities characterized by the same disorder strength and correlation length, i.e. the same 
localization length, but with different microscopic realization of the disorder (we assume 
that the pump profile -F(x) does not change from one disorder realization to the other). In 
the Anderson localized regime, CF states located deep enough inside the cavity are long- 
lived and therefore do not differ much from localized eigenmodes of closed systems. We thus 
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assume that they have the same statistical properties as closed cavity modes. Following 
Ref. [S5], we decompose d = 1 localized wavefunctions as 



iJm{x) = exp[-|x - Xm\/2^] 0„ (oj) / a/ 2^n(x) , (Sll) 

in terms of an exponential envelope and a real pseudo-random function (f)m{x) which remains 
correlated on the scale of the elastic mean free path i, {4>m{x)(l)n{x')) = 6mnGxp[—\x—x'\/2C^], 
where we used the fact that i = ^ for d = 1. 

In the Anderson localized regime, the correlation functions we need are 

{n\x)ijU^,k,)M^,k,)) = 5„„e-l^-^'"l/72e, (S12a) 

{n\x)n\x')^IJmix, k,)^m{x\ k^)Mx\ k^))m^n = e-(l--'"l + l-'-'"l+l--"l + l-'-"l)/2« 

xe-l^-^'l/«/4e^ (S12b) 
between eigenfunctions ipn{x, A;^) of a closed, Anderson localized system with unbroken time- 



reversal symmetry. Here, Xm is the center of mass of iprn{x). In Eq. (S12b) we neglect the 
hybridization contribution coming from even and odd linear combination of pairs of quasi- 
degenerate localized states [S6] but will comment on their possible effect below. In addition, 
localized modes have 

(|^^(x,A;^)|>„(x',A;^)|2) ^ {\i;^{x, k^)\^) x {\ij^{x',k^)\^) , (S13) 

when they are centered away from one another. 



From Eq. (S12a), we first conclude that only diagonal matrix elements Tmm have finite 
average, 

(T^n) = 5mn ^ dx exp [- |a; - | /^J 72^ . (S14) 
Unless one specifies F{x) and h{x), this is all one can say. 



We next evaluate var7^„. From Eq. (S12b) one has 



^ , , 1 w F'ix) 1 F{x') 1 

varTmn =11 dxdx 




c Jc 



1 + h{x) 'n?{x) 1 + h{x') rfiix') 



y. ^-{\x-X„i\ + \x'-X,n\ + \x-X„\ + \x'-Xn\)/2^ ^~\x-x'\/S, ^^^2 (S15) 

The integral vanishes, unless \xm — Xn\ < ^, and we substitute x„ — )■ Xm in Eq. (S15). The 



exponential exp[— |x — x'\/^] in the integrand in Eq. (S15) only leads to a prefactor of order, 
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but smaller than one, since the other exponentials already enforce that x and x' are within 
a localization length away from one another. We obtain the estimate 



(T 

\ I mm I 



(S16) 



for overlapping modes m and n. Note that %nm > always holds. 
We next give an upper bound for the average of 



A T2 



where corrections to the single-pole approximation are justified when 
in the main text]. We first note that Tmn = Tnm and thus 



(S17) 

< 1 [see Eq. (7) 



^ ^ -^nnTrnoriTrt 
nj^mo 



nmo 



I Ann I 5^r„ 

nj^mo 1 



2 

mgn ' 



(S18) 



We calculate 

/ ^ 'ri 



•2 

mgn 



dx 



dx J dx 

F{x) 
1 + h{x) 



, F{x) 



Fix'] 



1 + h{x) 1 + h{x') 
n2(a;) ' 



— ^rno {x)lpmo {x')lpn{x)tpn{x') 



(S19) 



where we used the completeness relation '^n'^"-{x)ilJn{x') = 6{x — x')/n'^{x). We therefore 
obtain 



^ max„^mo|A.,; 



dx 



^momo I 



F(x) 



/;,(x) 



n^(x) 



dx 



Ffx) 



n 2^ 



1 + /i(x; 



(S20) 



In the Anderson localized regime, the first factor on the right-hand side is exponentially 
small, max„^mJA„„|/|A„o^o = h/A^r^ ~ exp[-L/^], because (i) mo has kmo - and an 
exponentially large lifetime ~ {C,/c) exp[L/^], and (ii) the modes n ^ mo that overlap 



with ipmo differ in frequency by at least \u„ 



mo 



> A^/h, with the CF mode spacing 



hc/k'lj^ ^^'^ within one localization volume. The term between parenthesis, when averaged, is 
an algebraic function of /c^^, for any physically relevant choice of F{x), with the form of ipm 



given in Eq. (Sll). This is most easily seen, under the assumption that the pump profile 
F{x) = Fq is constant over the volume of the sample, which requires in particular the active 
medium to be homogeneous on the scale of the mode wavelength. When this is the case, we 



21 



obtain 



dx 



F{x) 



n 2 



1 + h{x) 



mo 



X) 



v?{x) 



SO that 



dx 



Fix] 



_l + h{x)_ 



Tl^ix) 



dx- 



P2 
-^0 



Fix) 



dx 
dx 



1 + 
1 



1 + /i(x 



mo 



n^i^x) 



(S21) 



< 1 



1 + /i(x) 



(S22) 



We expect this bound to still hold for smoothly fluctuating F{x). We conclude that is 
bound from above as < exp[— L/,^], up to algebraic corrections. 

This exponentially small upper bound on the average value does not rule out that S 
is large for some disorder realizations. We therefore investigate fluctuations of S. We obtain 
that they are still given by algebraic functions of k^^ times exp[— 2L/^]. The exponentially 
small factor comes from the pref actor max„^mo|A„„|/|Amomol- It can only be compensated 
if Tmomo — as exp[— L/,^] or faster. While this is in principle possible (e.g. for CF modes 
localized outside of the support of -F(x)), it never is the case for CF modes that contribute 
to lasing modes. We conclude that |S| ~ exp[— L/,^] ^ 1 for all disorder realizations in the 
Anderson localized regime. 



In the delocalized regime, Eq. (S20) still holds, and the term between parenthesis there 



becomes an algebraic function f{kf^i, fc^L) of k^L and k^£ [Eq. (S22) holds regardless of the 

hj /S.L'TmotJ, now depends 



lA 



momo I 



strength of the disorder]. The prefactor max„^mo|Ar 
on the level spacing = hc/27ik^L'^ and the mode lifetime Tmofi- The latter is no longer 
exponentially small. 

We close this paragraph with three comments. First, we recall that the infi- 
nite order nonlinearity arising from h{x) generates higher order correlators such as 
(^m(x, A;^)?/'„(x, kf,)[^^ T{k^) Y.p,p' <«i^Vp(x, ky)il)'pip^, K)]^), which the above treatment 
neglects, and which can be expected to be large for lasing modes. 

Second, Ref. [S6| points out that pairs (V'm,V'n) of quasi-degenerate localized wavefunc- 
tions exist that are spatially separated by more than ^. Because of their quasi-degeneracy, 
even and odd superpositions il)± = {ipm ± '4'n)/V^ are quasi-eigenmodes. They give 



rise to (i) a slower decay of the correlator in Eq. (S12b) at small frequency difference 



c{kmn — kn^) and (ii) a negative dip in the correlation function of Eq. (S12b) a distance 
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Lm — 2^1n[A^//j,(c/c^ju — cknn)]- If one such mode lases, the former effect is somehow rel- 
evant in that it enhances the ratio var7^^/(7Jnm)^ by a prefactor of order, but larger than 
one. The emergence of such Mott states, if they are lasing, can thus lead to an earlier onset 
of modal interactions. As for (ii) we believe that this effect is not relevant when F{x) is 
extended over more than Lm, because ipj^ and ^_ have equal but opposite correlations at 
X — Q and X ~ Lm- This mechanism for hybridization may however turn on modal inter- 
actions in Anderson localized random lasers when F{x) is localized over a region smaller 
than Lm- In particular it may potentially lead to lasing modes being spatially extended 
well outside the pump region. However, for such quasi-degenerate states, the static inver- 
sion approximation made to derive the SALT equations may break down, so that further 
investigations are necessary to determine the relevance of Mott states in random lasing in 
the Anderson localized regime. 

Third, we finally stress that from their very nature, all above estimates are parametric 
in essence, in that they can be multiplied by numerical prefactors of order one, depending 
on the exact form of F{x). 
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